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In this paper we investigate the influence of a constant external magnetic fleld on the flnite- 
temperature phase structure and the chiral properties of a simplifled lattice model for QCD. We 
assume an SU{2) gauge symmetry and employ dynamical staggered fermions of identical mass 
without rooting, corresponding to Nf = 4 flavors of identical electric charge. For flxed mass (given in 
lattice units) the critical temperature is seen to rise with the magnetic fleld strength. For three fixed 
/3-values, selected such that we stay (i) within the chirally broken phase, (ii) within the transition 
region or (iii) within the chirally restored phase, we study the approach to the chiral limit for various 
values of the magnetic field. Within the chirally broken (confinement) phase the chiral condensate 
is found to increase monotonically with a growing magnetic field strength. In the chiral limit the 
increase starts linear in agreement with a chiral model studied by Shushpanov and Smilga. Within 
the chirally restored (deconfinement) phase the chiral condensate tends to zero in the chiral limit, 
irrespective of the strength of the magnetic field. 
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I. INTRODUCTION 

It is commonly believed that hadronic matter at 
high temperatm'e undergoes a phase transition into 
another phase, traditionally called the "quark-gluon 
plasma" . Actually, ab initio numerical lattice simula- 
tions have shown that the high-temperature behavior 
of QCD at low baryon number density is governed 
by two interrelated crossover phenomena, namely the 
transition from a low temperature, confined regime to 
a high temperature deconfincd regime and the tran- 
sition from a low temperature regime with sponta- 
neously broken chiral symmetry to a high temperature 
regime with restored chiral symmetry. For a recent re- 
view see |l|]. The behavior at large baryon densities, 
on the other hand, is only predicted by models, be- 
cause in this case the lattice action is not real- valued 
in real QCD, which makes ab initio numerical sim- 
ulations impossible. There are special model theo- 
ries, however, where there is no sign problem at finite 
baryon density. This is the case, e.g., for SU{2) gauge 
theory with dynamical fermions. 

A very interesting question is how a strong external 
magnetic field modifies the properties of strong in- 
teractions at high temperature. Cosmological models 
suggest that very strong magnetic fields {sfeB ~ 1 — 2 
GeV) have been produced at the electroweak phase 
transition JZ]. Strong external magnetic fields cannot 
be ignored in noncentral heavy ion collisions at high 
energy, which are produced by the electric currents 
of the throughgoing spectator nucleons. Estimates of 
these fields range from \feB ^100 MeV for collisions 
at RHIC to \feB ~ 500 MeV at the LHC |3|, |Jj. 

One spectacular consequence in noncentral heavy 



ion collisions is the so called chiral magnetic effect, 
which was proposed in Ref. Q- In the presence of 
a winding number transition in the strongly interact- 
ing (color SU(3)) gauge field, the magnetic field can 
induce a parallel electric current. This leads to event- 
by-event fluctuations of the electric charge asymmetry 
of emitted hadrons with respect to the reaction plane. 
Such charge fiuctuations have, in fact, been seen in 
the STAR experiment at RHIC [1 Q and in the AL- 
ICE experiment at the LHC [7| . Whether these charge 
fluctuations actually come from the mechanism men- 
tioned above, is still under debate [8-10]. 

In the confinement region, the presence of an exter- 
nal constant magnetic field is expected to enhance the 
chiral symmetry breaking, i.e. to lead to an increase 
of the chiral condensate. This has been predicted in 
the Nambu-Jona-Lasinio model, as well as in the chi- 
ral model. Quantitatively, the growth of the chiral 
condensate with the magnetic field strength differs 
between the two approaches: in the NJL model the 
increase is predicted to be quadratic |lll - [l3| . while in 
the chiral model it is predicted to be linear in the chi- 
ral limit, for not too strong magnetic fields |1J|. The 
latter result is due to the infrared divergences occur- 
ring order by order in the magnetic field in the limit 
when the pion mass goes to zero. In the full solution 
(la ] the divergence disappears, because the magnetic 
field itself acts as an infrared regulator. 

According to several model calculations for the fi- 
nite temperature transition, the presence of an exter- 
nal magnetic field leads to an increase of the transition 
temperature [1^ JJ7| . The effect of an external electric 
field has been studied within holographic studies only 
|17| . In a two-phase treatment of the chiral model 



[isl [Tg l the temperature of the chiral phase transi- 
tion is seen to decrease with increasing magnetic field, 
while the originally first order transition (suggested by 
a bag model) ends with vanishing latent heat, becom- 
ing a crossover. The possibility of a different effect 
of a magnetic field on the chiral symmetry restoring 
transition on one hand (strongly increasing transition 
temperature) and on the deconfining transition on the 
other (moderately increasing transition temperature) 
has been pointed out in Ref. [20]. In the framework 
of a quark-meson model coupled to the Polyakov loop, 
the common crossover without magnetic field is seen 
to split into separate phase transitions. 

In the chiral limit, above the transition one expects 
the chiral order parameter to vanish. In Ref. [ij| it 
is shown with the help of Dyson- Schwinger equations 
for zero temperature in an external constant magnetic 
field (which is parametrically much larger than the 
other scales in the problem) that a quark mass propor- 
tional to •\/|ei?| is dynamically generated. This leads 
to a finite chiral condensate proportional to jeSp^. 
If these considerations can be applied to the high tem- 
perature phase, too, even in the chiral limit the phase 
transition would disappear in strong enough magnetic 
fields. 

Recently several lattice simulations have been per- 
formed in gauge theories with fermions coupled to a 
constant external magnetic field. ^ The inclusion of the 
magnetic field does not lead to a sign problem for dy- 
namical fermions. The pioneering work was performed 
in Refs. [22 - l25J |. There the simulations were done 
in quenched SU{2) gauge theory coupled to overlap 
fermions. The chiral condensate was found to increase 
linearly with the magnetic field strength in the chiral 
limit even in very strong magnetic fields \/eB < 3 
GeV, both at zero temperature and at T = 0.82Tc 
slightly below the critical temperature Tc- A later in- 
vestigation, still within quenched QCD, tends to favor 
an increase of the chiral condensate with the magnetic 
field strength with an effective power 1.6(2) [26]. Sim- 
ulations in full QCD with two fiavors using staggered 
fermions have been reported in Refs. [27|,|2^. The au- 
thors of these works find, in agreement with the ear- 
lier quenched simulations, that the chiral condensate 
grows with the magnetic field strength, but quadrati- 
cally for not too large magnetic fields. They attribute 
this to their use of a finite quark mass, corresponding 
to a pion mass tti^ « 200 Mev. They further find that 
the chiral condensate increases with the magnetic field 
in the whole range of finite temperatures, thus even in 
the transition region and beyond [27|. This, in turn. 



leads to an increase of the finite temperature tran- 
sition temperature with the magnetic field strength. 
An opposite conclusion was very recently presented 
in [29|, namely that the transition temperature de- 
creases with an increasing magnetic field strength. 
This comes about because the chiral condensate, al- 
though increasing in the confined region away from 
the transition, actually decreases with the magnetic 
field strength in the transition region. In this inves- 
tigation the authors have performed the simulations 
using an improved action with 2-1-1 fiavors of dynami- 
cal staggered fermions. They suggest that their results 
(differing from others) are due to the fact that in their 
calculation the pion mass is lower, and the taste split- 
ting considerably smaller. Both calculations use the 
"fourth root trick" for a single fiavor in order to rep- 
resent the right number of (non-degenerate) flavors. 

In this article we investigate the case of SU{2) gauge 
theory with dynamical fermions at finite temperature. 
The influence of the constant magnetic field should to 
a large extent depend on the chiral properties, as pro- 
posed in Ref. [1J|. These are quite similar to those 
of QCD, even if the color group is different. Further- 
more, in the SU{2) theory one can extend the investi- 
gation to finite chemical potential, as well as study in 
detail the topological excitations, which should be re- 
sponsible for the chiral magnetic effect. In this paper 
we still confine our interest to the response of the chi- 
ral condensate and the finite temperature transition 
to the magnetic field. 

We simulate the theory on the lattice, using stag- 
gered fermions, without invoking the root of the 
fermion determinant. In the continuum limit this 
leads to a four-fiavor theory with equal electric charge, 
but we avoid the bias due to rooting. We can de- 
termine the influence of the unquenching by compar- 
ing with the simulations in quenched SU{2) theory in 
Refs. [2a - t25| . We can, of course, also test the validity 
for this case of the models mentioned above. 

In Section In] we define the action and the order pa- 
rameters. In Section IlIII we describe the setup of our 
simulations and how we determine the scale by calcu- 
lating the heavy quark potential and the pion mass at 
zero temperature. Section IIVI presents the results at 
finite temperature and with magnetic field for the chi- 
ral condensate, the Polyakov loop and the mean values 
of plaquettes differently oriented with respect to the 
magnetic field. Section |V] contains the conclusions. 



II. SPECIFICATION OF THE ACTION AND 
DEFINITION OF THE ORDER PARAMETERS 



^ We emphasize that this task is conceptionally different to 
those where an external co/our-magnetic field acts addition- 
ally to the fluctuating colour gauge field, see e.g. |2f| for a 
recent work. 



To describe the finite temperature theory we intro- 
duce a lattice of size Nr x N^ . The sites are enumer- 
ated by n = (ni, ri2, "-3, "•4), where the rit are integers. 
Hi = 1,2,..., Ncr for i = 1, 2,3 and 714 = 1, 2, ... , Nr- 
The fourth direction is the (inverse) temperature or 



(Euclidean) time direction. The lattice spacing is de- 
noted by a. On the links n — >■ n + /i we define group 
elements U^{n) € SU{2), where n = 1,2,3,4. For 
the gauge action we adopt the usual Wilson plaquette 
action 



/3 



^G = ^EEt^-(i-t^A..w) 



n /i<i/ 



where U^u{n) is the ^v plaquette matrix attached to 
the site n: 

U^,,{n) = U^{n)U,{n + fi)Ul{n + i>)Ul{n). (2) 

We want to introduce an external constant magnetic 
field to interact with the fermions. We therefore in- 
troduce electromagnetic potentials in the fermion ac- 
tion by new, commuting group elements on the links, 
namely Vfj,{n) — e'^f*^") S U{1), with compact link 
angles < 0^{n) < 2n. We further introduce stag- 
gered fermions as Grassmann variables lij^n) and ip{n), 
which are color vectors in the fundamental represen- 
tation of SU{2). The fermionic part of the action 
becomes 

Sf = a3^V^(p)[i^(p,(7)+ma5p,,]V(g), (3) 
p,q 



- V:ip-ii)Ulip-^,)S^^ 



P-M:9J 



(4) 



The arguments p, q are integer four vectors denoting 
sites on the lattice and /^^(p) are the usual staggered 
sign factors, 

771 (n) = 1, (5) 

ry^n)- (-1)^-1 "^ /i = 2,3,4. (6) 

In Ref. [30| the authors propose a construction to in- 
clude a constant magnetic field both in the continuum 
and on the lattice, with periodic boundary conditions 
in the spatial directions. In the continuum, to have a 
constant magnetic field B — (0, 0, B) pointing in the 
z-direction they define the vector potentials as 

Af,{x,y,z,t)^ —{xSf,^2~ySf,,i)- (7) 

Because of the periodic boundary conditions, delta 
functions are needed at the boundary. Then B = 
rot A is constant except on the boundary, where we 
obtain a large magnetic field, so that the average of 
the magnetic field is zero 



1 



LxLy 



dxdy Bz 



0. 



(8) 



' (a;, y) — plane 

The translation of this construction to the lattice is a 
good choice, because the lattice sees the delta function 
only modulo 27r. A plaquette angle can be defined by 



V(n)=A^6i^-A^6i^ 



(9) 



where A^f{n) = {f{n + u) — f{n))/a is the lattice 
forward derivative acting on 9y{n) ^ A^(n). Then 
the electromagnetic plaquette part can be split as [3l| 



[0]f^,yin) +2Tm^i,(n) 



(10) 



(\\ with n^^{n) = —2,-1,0,1,2. Plaquettes with 71^,^ ^ 
are called Dirac plaquettes, and the reduced pla- 
quette angle [0]^^ G [0,27r[ corresponds to the (gauge- 
invariant) electromagnetic flux through the plaquette. 
A constant magnetic background field in the z- 
direction penetrating all the (a;, y) -planes of finite size 
No- X N^ with a constant magnetic fiux (j) per each 
(x, y)-plaquette is realized by the following choice: 



Vi{n) 



-i4>n2/2 



(ni = l,2,. 



■ ,N„ 



1), (11) 



^2(71) = e^*"i/2 (n2 = l,2,...,7V,-l), (12) 

Fi(iV,,n2,n3,n4) = e-*'A(^'+i)"^/2^ (13) 
F2(ni,7V„n3,n4) = e*^(^-+i)"i/2_ (^4) 

In order to have a constant flux also on the boundary 
it has to be quantized as follows {q denoting the same 
electric charge of all fermion flavors). 



a^qB 



2TTNb 



mez. 



(15) 



Thus, on the lattice there is always a minimal non- 
vanishing flux. Because the angle is periodic, there is 
also a maximum flux. In fact one has to restrict oneself 
to (/) < TT, or Nb < N^/2. At flnite temperature this 



means that -^4?^ is restricted to the region 



-Nr 



yptix-L < yjL- < ^Nr. 



N, 



T 



(16) 



In the case of color group SU{2) there is a larger chiral 
symmetry than in color SU{N) with N > 2, because 
the fundamental representation is equivalent to the 
conjugate one. In the continuum, for B — m — 
there is a U{2Nf) chiral symmetry, which is broken 
down to Sp{2Nf) when the mass is different from zero. 
If the symmetry is spontaneously broken, there are 
Nf{2Nf — 1) — 1 Goldstone bosons, where the —1 (one 
less Goldstone boson) is due to the axial anomaly. 

On the lattice, for SU{N), N > 2 there is a global 
C/(l) (8) f/(l) chiral symmetry for i? = m = 0. The 
nonsinglet axial symmetry is broken when m 7^ 0. 
Lattice calculations show that it is also spontaneously 
broken, with one Goldstone boson. In color SU{2) 
ior B = m = there is instead a global U{2) chiral 
symmetry on the lattice. For m ^ 0, B = it is 
broken down to f7(l). In the dynamical theory it is 
spontaneously broken giving rise to three Goldstone 
bosons. One of them is an electrically neutral meson, 
while the other two are a baryon and an antibaryon, 
which are electrically charged. Notice that in color 
SU{2) the baryons are bosons. When B ^ there is 



only one Goldstone meson. The baryons get masses 
proportional to ^/qB. The situation, as far as the 
Goldstone bosons are concerned, is thus quite similar 
to two-flavor QCD. 

In order to study the influence of the external static 
magnetic field on the chiral condensate and on the 
phase structure we measure the following (approxi- 
mate) order parameters. The chiral condensate, which 
is an exact order parameter in the limit of vanishing 
quark mass, is given by 

119 

NrNl 4 d{ma) ^^ ' 



1 1 



NrNl 4 
where the partition function is 

Z= I Y[{dtP{n)dtP{n)dU^,{n))e 

71, fl 



- <TT{D + ma)-' >, (17) 



— Sg — Sf 



(18) 



The factor 1/4 is inserted because we define < ipip > 
per flavor, and our theory has 4 flavors. 

In order to locate the phase transition we use the 
disconnected part of the susceptibility (later on called 
"chiral susceptibility" for simplicity). 



1 1 



d^ 



log(Z) 



Xdu 



NrN^ 4 {d{ma))^ 

Xconn t Xdisc ; 

^l«(Tr(i. + ^a)-T> 



(19) 



(20) 



It is important to notice that these are bare quantities, 
which should be renormalized when comparing with 
continuum expectation values. 

We further measure the average value of the order 
parameter for confinement, the Polyakov loop 



<^>= Nl ^ 2 

n 1,712,^3 



(21) 



< Tr ]^ Ui{ni,n2,nz,ni) > 

\n4 = l J 



and the corresponding susceptibility 

XL = N^{< L^>-<L >2). (22) 

Intuitively it is clear that a constant magnetic field 
oriented into one of the three space directions has to 
violate the 3d isotropy of the system. The isotropy 
violation can only be caused by the coupling of the 
magnetic field to the fermionic part of the action and 
should become visible within the effective gauge action 
after the fermion degrees have been integrated out. 



In order to demonstrate this anisotropy one can eas- 
ily compute the average non-Abelian plaquette values 
separately for the different space-time planes (p,, v) 
(an averaging over the lattice is implied) 



Pt^u =< \ Re Tr U^, > , 
as a function of the magnetic flcld strength. 



(23) 



III. SETUP OF THE SIMULATIONS 

In order to simulate the model with dynamical stag- 
gered fermions (Nf = 4) we employed the standard 
Hybrid Monte Carlo algorithm. We chose the number 
of integration steps Hint and their "time" length St 
such that the length of a trajectory was nintSr = 1 
while the acceptance rate was larger than 0.8. For the 
finite-temperature measurements we used a 16"^ x 6 
lattice. We measured the chiral condensate on ev- 
ery fifth configuration, the Polyakov loop and the pla- 
quette variables on every configuration. The chiral 
condensate was calculated with the random source 
method. Thereby we used 100 random sources per 
configuration. The integrated autocorrelation time 
was taken into account in all our error estimates. 
In general, apart from simulating very near to the 
transition temperature, the integrated autocorrelation 
times of all observables were estimated mostly well 
below 20. The number of configurations (trajectories) 
generated within a run varied between 1800 and 5000. 
In general, 300 configurations were discarded for ini- 
tial thermalization. 

We also made a scan of the susceptibilities Xdisc and 
XL near the transitions for different magnetic fields. 
For this aim we generated longer runs of length be- 
tween 10'* and 2 • 10^ trajectories, while we measured 
the chiral condensate and the Polyakov loop with the 
usual frequency. The autocorrelation times were esti- 
mated to 0(25) near the transition points. 

Some zcro-tempcrature simulations on a lattice of 
size 16"^ X 32 without magnetic field were performed 
additionally in order to estimate the lattice spacing, 
the pion mass, and the critical temperature in physical 
units. We are, of course, aware of the fact that we are 
considering a fictitious world of two-color QCD with 
four flavors of identical electric charge q. Nevertheless, 
a scale determination provides a rough estimate how 
near we are to the chiral limit, and how large the 
explored magnetic fleld strengths are. The number 
of trajectories per /3-value used in the T = case 
was 2000. Measurements were performed after every 
fourth trajectory. 

At T = we first calculated the pion propagator 
C7r(i) for the three /^-values 1.7, 1.8 and 1.9, fitting it 
to the usual form (Nr — 32) 



a(f)=Ce-^* + Ce^(*-^^). 



(24) 



The fit parameters in lattice units can be found in 
Table HI Our results are compatible with earlier ones 
obtained for larger quark masses on smaller lattices 

Moreover, we computed the heavy quark poten- 
tial from HYP-smeared [3j] Wilson loops (with HYP- 
smearing in the version of R ef. l34l| ) and for compari- 
son also with APE-smearing [3^ at the three /3- values 
1.8 , 1.9 and 2.1, by fitting it to the form 



with 



V{R) = Vs{R) + C{--Gl{R)) 



Vs{R) =A-B/R + aR, 



(26) 



where Gl{R) is the free gluon propagator on the lat- 
tice. This correction is only important at short dis- 
tances. For details of this procedure see W]. The 
fit results obtained from the HYP-smeared data are 
presented in Table HI too. 

In order to estimate the lattice spacing in physical 
units we have determined the Sommer scale parameter 
Ro [37] related to the static force Fs{R) = dVs/dR 
through the condition 



Fs{Ro)Rl = l.65. 



(27) 



The values Ro/a found at ma = 0.01 for the three 
/3- values (1.8, 1.9, 2.1) are given in Table HT] To get 
the scale in physical units, we adopted the physical 
value ro = 0.468(4) fm fixed for QCD [38]. Our mea- 
surements of the scale a at the three /?- values are well 
fitted by the two-loop formula for the /3-function for 
the SU{2) case and four flavors 



a(/3) 



1 ,4/3, 



213, 



Al' (3 



(^)23fi exp(-f) 



A 



24 



I- 

19 

384^ 



(28) 



with Al = 0.00660(12) fni"^ and xL/ = 1-64 (see 
Fig. [T]). Thus, we used this formula for an extrapola- 
tion down to f3 = 1.7, where the lattice spacing a is 
too large to allow a safe determination of Rq. In this 
way we could roughly estimate the pion mass and the 
magnetic field strength also at this value of the bare 
coupling and ma = 0.01. In Table lU we have collected 
all our results for setting the scale. 

Let us mention, that for finite temperature with 
iV^- = 6 and ma = 0.01 (for vanishing magnetic field 
strength) the critical value f3c will be seen close to 
/? = 1.8. Thus, we are able to estimate also the critical 
temperature in physical units. We find Tc ~ 193(6) 
MeV. 



IV. RESULTS AT FINITE TEMPERATURE 

The finite temperature simulations are performed 
on lattices of size 16^ x 6. In Fig. [5] the bare chiral 



condensate is shown as a function of the fiux num- 
ber Nb for various values of /3 in the region of the 
transition value /3 = /^^ = 1.8. Because the lattice 
is finite the chiral condensate is periodic in the mag- 
netic flux (j). We see that saturation effects set in at 
Nb « 60. This corresponds to Nb « N^/A or ss 7r/2, 
which is half of the maximally achievable fiux, see for- 
mula (|16p . Our investigation is mostly restricted to 
values of Nj, < 50. The corresponding maximal val- 
ues of \/qB are given in Table |lll They are in general 



(25) greater than the physical region of interest. 



GeV. In Fig. [3] we plot the bare chiral condensate as 
a function of /3 for a set of numbers of flux quanta 
for the bare quark mass ma = 0.01 and 0.1, respec- 
tively. For the smaller quark mass we see quite clearly 
a transition for all values of the flux quanta Nb under 
consideration. It is important to notice that the bare 
chiral condensate increases with increasing magnetic 
field for fixed /3, for all values of /3 in the transition re- 
gion. This indicates that the chiral transition moves to 
higher temperatures as the magnetic field is increas- 
ing. This tendency is in agreement with the results 
in [23, [2I] but opposite to the tendency seen in [20] 
where the bare chiral condensate decreases with the 
flux (f> in the transition region, leading to a decrease 
of the transition temperature. As shown in Section 
III, without magnetic field, for am, = 0.01 and close 
to Tc we have reached a ratio ttIt^/Tc ~ 1.7, which is 
similar to the ratio in [23, [2^ , but higher than that 
in [20]. We also have a different gauge group. It is 
clear that our observation does not represent a direct 
contradiction to the results in |29| . 

For the higher quark mass ma = 0.1 the transition 
in the chiral condensate turns out very smoothly. 

In Fig. |4] the expectation value of the Polyakov loop 
is shown vs. /3 for the same two values of the bare 
quark mass. This is an indicator for the deconfine- 
ment transition. For the lower quark mass the decon- 
finement transition and the chiral transition are con- 
sistent with happening at the same temperature. Fur- 
thermore, the transition temperature increases with 
the quark mass. At the high quark mass there seems 
to be only a weak effect of the magnetic field on the 
deconfinement temperature. 

To study the change of the transition temperature 
at the lower quark mass more quantitatively we have 
calculated the susceptibilities. In Fig. [5] we show the 
chiral susceptibility and the Polyakov loop suscepti- 
bility for the low quark mass ma = 0.01. It is clearly 
seen in the left figure that the chiral transition in- 
deed moves to higher temperatures as the magnetic 
field becomes stronger. In the right figure we show 
the same effect for the Polyakov loop susceptibility. 
In fact, the maxima of the two susceptibilities are at 
the same value for given magnetic field. There is no 
sign of a splitting between the chiral and the decon- 
finement transition. The height of the peaks increases 
with the field strength. However, only a finite size 
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E — arriT^ 
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2 
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1.7 
1.8 
1.9 
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1.48(5) 

1.01(3) 

0.311(11) 


0.265(1) 
0.285(1) 
0.296(2) 


0.047 
0.023 
0.039 


0.265(38) 
0.265(9) 
0.250(9) 


0.370(42) 
0.300(11) 
0.265(13) 


0.169(8) 
0.0725(19) 
0.0170(16) 


0.872 
0.681 
1.93 



TABLE I: Fit parameters in lattice units for the pion correlator C,r(i) ace. to Eq. (|24p and for the static potential Vs{R) 
ace. to Eq. (|26|l . The fit range for the pion correlator starts at lattice distance tmin- 
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Ro/a 


a[fm] 


m^MeV] 


V^^^^solGeV] 


1.7 


1.89(4) 


0.248(4) 


210(4) 


0.881(15) 


1.8 


2.75(8) 


0.170(5) 


330(10) 


1.29(4) 


1.9 


4.32(6) 


0.108(2) 


537(9) 


2.02(3) 


2.1 


9.03(43) 


0.052(2) 


- 


4.22(21) 



TABLE II: Results in physical units for the Sommer scale Ro, the lattice spacing a, the pion mass m.^-, and the quantity 
y/qB characterizing the magnetic field strength for Nt — 50 flux units for various /?- values considered in Section HVl The 
values for /? — 1.7 were estimated by extrapolating with the two-loop beta-function. 



scaling analysis could show if the transition remains a 
crossover or becomes a real phase transition. 

In order to study the dependence of the chiral con- 
densate on the magnetic field strength in the chiral 
limit, we have chosen to investigate in more detail the 
behavior of the chiral condensate as a function of the 
quark mass and the magnetic field for three fixed val- 
ues of /?. Because we keep (3 fixed as we vary the 
quark mass and the magnetic field, the lattice spac- 
ing a is also fixed and we have eliminated lattice ef- 
fects coming from the variation of a. The three val- 
ues are /? — 1.70 (which is clearly in the confined 
phase), /? = 1.90 (which is in the transition region), 
and /3 — 2.10 (which is in the deconfined phase). Let 
us first discuss the results in the confined region. 

In the left panel of Fig. [6] we show the dependence of 
the bare chiral condensate on the quark mass for vari- 
ous values of the magnetic flux. To obtain the results 
relevant to continuum physics, one has to subtract an 
additive divergence for finite quark mass, as well as 
do a multiplicative renormalization, which is needed 
also at zero mass. In the right panel of Fig. [S] we show 
the difference between the bare chiral condensate for 
finite fluxes subtracted by the same quantity at zero 
flux. This eliminates the main part of the additive 
divergence. In the left panel we have also included 
points at quark mass zero, where there are no addi- 
tive divergencies. The values at zero quark mass are 
obtained by a chiral extrapolation. We perform this 
extrapolation in two ways. Because we are not very 
far from the transition, we suppose that we can use 
the formula for the reduced three dimensional model 
[39j . See also [40l - |42| . In that case we may use the 
ansatz a^ < ipip >= fi{ma), where 



with the non-analytic term coming from the Gold- 
stone bosons. Such a parametrization has been used 
in [38| in the context of the finite temperature transi- 
tion in QCD. As an alternative we also use the chiral 
extrapolation relevant for zero temperature, namely 



< f/'V-' >= /2("T-a), where 

/2(ma) = 6o + 6imalogTOa + b2ma. 



/i (ma) = ao 



a2ma 



(29) 



(30) 



The fits performed with the Eq. (f29|) are shown as 
dashed or dotted lines in Fig. [T] The fit parameters 
of both kinds of fits are summarized for /3 = 1.70 and 
various magnetic fiuxes in Table IIIII 

The value of the bare chiral condensate in lattice 
units within the chiral limit (m = 0) for various values 
of the flux are given by oq and bo, respectively. One 
can see that the statistical errors on the parameters 
are much smaller than the systematic errors coming 
from the fit formulas. We are, however, not interested 
in the absolute values of the bare chiral condensate at 
zero mass, but in the dependence of this value on the 
magnetic field strength. The flux dependent part does 
not need renormalization, and thus can be directly 
compared to continuum models. In Fig. [7] we show 
the bare chiral condensate as a function of the num- 
ber of magnetic flux quanta for different bare quark 
masses. The chirally extrapolated points are obtained 
with the fit functions in Eqs. ([29|) and (pO| . One can 
see that for finite quark mass the data at small mag- 
netic flux are consistent with a quadratic behavior, 
while the extrapolated values to ma = seem to start 
with a linear behavior irrespective of the specific chi- 
ral extrapolation /i or /2 used. This behavior is in 
agreement with the prediction of the chiral model of 
Ref . [IJ] . In Table IIVI we provide the fit parameters 
for the chiral condensate as a function of the magnetic 
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Nb 
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Xdof 



2 

10 
20 
30 
40 
50 


0.080(1) 
0.083(1) 
0.098(2) 
0.1218(8) 
0.1375(4) 
0.1492(6) 
0.1586(7) 


0.36(1) 
0.34(1) 
0.28(1) 
0.203(3) 
0.177(1) 
0.162(2) 
0.149(2) 


-0.19(2) 
-0.16(2) 
-0.09(3) 










0.33 
0.20 
0.66 
0.73 
0.20 
0.35 
0.23 


0.098(1) 
0.100(1) 
0.1118(5) 
0.1314(2) 
0.1459(3) 
0.1568(2) 
0.1657(5) 


-0.42(1) 
-0.42(1) 
-0.28(1) 
-0.238(1) 
-0.208(2) 
-0.1901(9) 
-0.174(2) 


-0.20(2) 
-0.16(3) 
0.08(1) 










0.73 
0.30 
0.10 
0.80 
0.18 
0.05 
0.18 



TABLE III: Chiral fit parameters for the fits fi{ma) (Eq. (|29|) ) and f2{ma) (Eq. pO[) ) allowing to extrapolate to the 
chiral limit for /3 — 1.70 (chirally broken phase) and various magnetic field strengths. The fit curves obtained with /i 
are shown in the left panel of Fig. [6] 





giNt) 




Co 


Cl 


C2 


2 
Xdof 


CEl 
CE2 


0.080(1) 
0.095(2) 


0.027(2) 
0.021(2) 


-0.00015(5) 
-0.00015(3) 


1.18 
9.4 



TABLE IV: Fit parameters for the chiral condensate as a function of the magnetic flux A^i, in the chiral limit ace. to 
the polynomial ansatz Eq. (|3ip . CEl and CE2 correspond to the chiral fit extrapolations with fit functions /i and /2, 
respectively, as presented in Table IIIII Both the fit results are shown in Fig. [7] 



flux (j)^ Nb 



g{Ni) = cn{l + ciNb + C2NI) 



(31) 



in the chiral limit obtained with the chiral extrapo- 
lation /i (CEl) and /2 (CE2), respectively. We have 
not calculated the pion decay constant, and therefore 
we do not make a quantitative comparison with the 
chiral model at finite temperature and magnetic field 
studied in Ref. [l8l| . 

Now we come to the transition region and the tem- 
perature region above the transition. In Fig. [5] we 
show the mass dependence of the bare chiral conden- 
sate (left panel) and the subtracted chiral condensate 
(right panel) in the transition region (at /3 = 1.90) for 
three values of the magnetic flux. One can see that 
for finite flux, as well as for zero flux, the bare and 
subtracted chiral condensates are consistent with ex- 
trapolating to zero in the chiral limit. For the highest 
flux, iV(, = 50 one can clearly discern two regions of 
behavior. For aTO^0.04 the chiral condensate seems 
to extrapolate to a finite value, but for am^0.04 it ac- 
tually extrapolates to zero. This can be understood, 
if one assumes that the transition for N^ = 50 at this 
value of P takes place for am « 0.04. In Fig. IH] we 
present the same quantities as above, but for /3 = 2.10. 
This is well inside the chirally restored phase. The 
chiral condensate extrapolates to zero for all values of 
the flux. Thus chiral symmetry is restored for all val- 
ues of the flux that we have investigated. Within our 
precision there is no evidence for a quark mass which 
is spontaneously created by the magnetic field as sug- 
gested in a selfconsistent calculation in [14] using the 
Dyson-Schwinger equations. 



Finally, we study the indirect influence of the mag- 
netic fleld on the gluonic part of the action. For this 
purpose we provide the plaquette energies 1 — P^,y in- 
dividually for the six differently oriented plaquettes, 
as given in Ea. (j23p for varying fermionic mass. We 
do this for the same values of (3 as above, namely 
/3 = 1.70, 1.90 and 2.10 and for two values of the flux 
corresponding to Nf, = (zero magnetic fleld) and 
Nh = 50 (strong magnetic fleld). 

In Fig. [TOlone can see that at /3 = 1.70 (within the 
conflned phase) all the plaquettes are degenerate for 
vanishing magnetic field. For a strong magnetic field 
pointing into the z direction, the xy and the zt pla- 
quettes become well separated from the others. This 
effect can be qualitatively understood by represent- 
ing the effective gauge action in terms of fermionic 
loops within a hopping parameter expansion for large 
fermion mass am ^ 1/k. While in the order «;■* of 
such an expansion only xy-plaquettes receive a cou- 
pling to the magnetic field strength, in the order k^ 
loops extending into three directions xyz or xyt may 
couple to it. The zi-plane is distinct, since only in the 
order k® and beyond there exist loops extending into 
all four directions which couple to the magnetic field 
in the effective action. 

In Fig.[TT]the same plots are shown for j3 — 1.90 (in 
the transition region). For vanishing magnetic field 
now the spacelike and the timelike plaquettes are dif- 
ferent from each other for sufficient small quark mass. 
This is the well-known temperature effect observed in 
the deconfined phase and providing a non-vanishing 
energy density (see e.g (43|). However, in the right 
panel one can see that, with a strong magnetic field 
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applied, the splitting of the xy and the zt plaque- 
ttes from the others is even stronger than the splitting 
between the other spacelike and timelike plaquettes. 
Deep in the deconfined phase these results become 
even more pronounced, see Fig. 1121 



V. CONCLUSIONS 

We have investigated two-color QCD at finite tem- 
perature in an external magnetic field using lattice 
simulations. We have in particular studied how the 
magnetic field infiuences the chiral properties of the 
theory. As the chiral properties of this theory are sim- 
ilar to those of QCD, our results should be relevant 
also for the physical case. 

We have found for all temperatures for fixed bare 
quark mass that the chiral condensate grows with the 
magnetic field. Hence, the temperature of the chi- 
ral phase transition grows with the strength of the 
magnetic field. This is confirmed by a measurement 
locating the peak of the chiral susceptibility. At the 
phase transition the pion mass is around 1.5 times the 
critical temperature. A similar result had been found 
by t23) ^M i^ QCD with about the same ratio of the 
pion mass to the critical temperature. The opposite 
conclusion was reached in [29| , claiming that the tran- 
sition temperature decreases with the magnetic field 
in QCD with the physical pion mass. It would be in- 
teresting to see if the result in our model is valid in 
the chiral limit, but we leave this to a future inves- 
tigation. Furthermore, we have shown, by measuring 
the susceptibility of the Polyakov loop that the de- 
confinement and the chiral transitions move together 
to the same temperature, if small or large magnetic 
fields are switched on. 

In this work, by extending our measurements to sev- 
eral values of the quark mass, we make a first investi- 



gation of the chiral limit of the theory at three values 
of the temperature, one in the confined region, one in 
the transition region and another deep in the decon- 
fined phase. In the chiral limit there is no additive 
renormalization necessary for the chiral condensate. 
In the confined region, we find that the chiral conden- 
sate in the chiral limit seems to grow linearly with the 
magnetic field strength. This is in agreement with the 
prediction of the chiral model at zero temperature. 

At the other two values of the temperature chosen, 
the chiral condensate extrapolates to zero for all val- 
ues of the magnetic field. Thus, in the chiral limit 
there is a real chiral phase transition, which does not 
disappear for strong magnetic fields. 

We have also investigated the influence of the mag- 
netic field on the gluonic part of the theory. We find an 
asymmetry of the non-Abelian plaquette values with 
respect to the magnetic field. These variables are re- 
lated to the gluonic part of the energy density of the 
theory. Thus we expect an influence of the magnetic 
field also on the equation of state. It would be inter- 
esting to conduct a study dedicated to this effect. 
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FIG. 1: The lattice spacing a vs. /3. The two-loop beta-function ace. to Eq. (|28p was used as fit-function, and a(/3 = f .7) 
was obtained through extrapolation. 




FIG. 2: The rising and saturation behavior of the chiral condensate as a function of the magnetic field in flux units for 
am = 0.01 and for various /3-values. The lattice size is 16^ x 6. The curves are to guide the eyes. 
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FIG. 3: The bare chiral condensate a^{'4>'4>) vs. inverse coupling /3 for various magnetic fluxes 4> (in flux units) and for 
two different bare masses ma = O.OI (left panel) and am = 0.1 (right panel). The lattice size is f6^ x 6. The curves are 
to guide the eyes. 
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FIG. 6: Mass dependence of the bare chiral condensate (left) and of the subtracted chiral condensate (right) for various 
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FIG. 9: Same as in Fig. [S]but witliin the deconfinemeiit piiase (/? = 2.1). 
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FIG. 11: Same as in Fig. [TOlbut for /? = 1.90 (transition region). 





FIG. 12: Same as in Fig. 1101 but for /3 = 2.10 (deconfinement phase) 



